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Abstract 

This paper is devoted to study the analytic models of anisotropic 
compact stars in f(T) gravity (where T is torsion scalar), with non¬ 
diagonal tetrad. By taking the anisotropic source inside the spheri¬ 
cally symmetric star, the equations of motions have been derived in 
the context of f(T) gravity. Krori and Barua metric which satisfies 
the physical requirement of a realistic star, has been applied to de¬ 
scribe the compact objects like strange stars. We use the power law 
form of f(T) model to determine explicit relations of matter variables. 
Further, we have found the anisotropic behavior, energy conditions, 
stability and surface redshift of stars. Using the masses and radii of 
4U1820-30, Her X-l, SAX J 1808-3658, we have determined the con¬ 
stants involved in metric components. Finally we discuss the graphical 
behavior of the analytic description of strange star candiddates. 
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1 Introduction 


The accelerated expansion of the Universe meg, and its verification by sev¬ 
eral observations E-0, motivate the theocratical physicists to propose such 
theories of gravity, which can extend General Relativity (GR). The extension 
should be in such a way that observation can be justified by the modifications 
instead of introducing dark energy. It is well known that GR deals with the 
description of gravity at local scale but not at global scale, for example GR 
fails to explain the the rotation of spiral galaxies |6]. Several theories of grav¬ 
ity have proposed to generalized the Einstein theory of gravity, but most of 
these theories are the geometric extension of GR which provide the essential 
ingredients to support the observations. Among all these modifications, one 
can consider the f(R) theories, in which Lagrangian can be expressed as a 
function of Ricci Scalar R [7, S]. 

Analogous to f(R) theory of gravity there is another theory so-called 
f(T ) theory, this theory is defined by the teleparallcl equivalence of gravity 
(TEGR). In this theory gravitational action can be described by the torsion 
scalar T instead of Lagrangian. In f(T) theory, Riemann-Cartan spacetime 
is defined with Weitzzenbbck connections instead of Levi-Civita connections, 
in this situation background spacetime admits non-zero torsion and zero 
curvature. This definition of spacetime was introduced by Einstein to define 
gravitation in terms of torsion and tetrad. In TEGR, tetrad field plays the 
role of dynamical field instead of metric field. Inspite of such differences 
TEGR and GR have many similarities, in other words every solution of GR 
is also solution of TEGR. But when TEGR is generalized to /(T) gravity 
theory, by assuming the Lagrangian as a general function of torsion scalar, 
the equivalence between these tow theories breaks down mum Therefore 
f(T) theory can be considered as a prominent candidate for the explanation 
of the accelerated expansion of universe in the absence of dark energy [T3] . 

In f(T) gravity, equations of motions are second order differential equa¬ 
tions as in GR, but in f(R) such equations are fourth order in metric for¬ 
mulism (second order in plantini approach). As compared to GR, dynamical 
equations in f(T) gravity, displays the additional degree of freedom which is 
related to the fact that equations of motions are invariant under the Lorentz 
transformation [15] . Recently [IB], it has been proposed that there are two 
types tetrad, bad and good. The diagonal and non-diapasonal tetrad are 
known as bad and good tetrad, respectively. The non-diagonal tetrad are 
used to properly parallelize the isotropic spherically symmetric spacetime. 
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No doubt f(T) theory is excellent theory for the explanation of cosmic ac¬ 
celeration and observations on large scale(clnstering of galaxies) [17]. Bnt 
we must remember that GR is in excellent agreement in with solar system 
test and pnlsar observation [IS]. In recent papers [19, .20], some solar system 
constraints have been proposed for the f(T) gravity in the diagonal tetrad 
case. Bnt in the present paper, we adopt the non-diagonal tetrad to drive 
the models of anisotropic compact stars in f(T), gravity by assuming the 
power law form of f(T) model. 

In theoretical astrophysics, f(T) theory was used to explore the effects of 
f(T) models in 3dimensions, so that f(T) version of BTZ black hole solutions 
has been formulated [21]. Later on [ 22 ], it has been proved that the violation 
of first of black hole thermodynamics in f(T) gravity, is due to the violation of 
Lorentz invariance. Recently [23], some static spherically symmetric solutions 
with charged source have been found in f(T) theory. After finding a large 
class of static perfect fluid solutions [23], the physical conditions for the 
existence of astrophysical stars in f(T) theory have been discussed [25]. 

Recently, Abbas and his collaborators have studied the models of anisotropic 
compact stars in GR, f(R), f{G ) and f(T) theories in diagonal tetrad case 
[26]-|29j. The main objective of this paper is to study the symmetric models 
of the anisotropic compact stars in the context of f(T) gravity using non¬ 
diagonal tetrad. We investigate the anisotropic the anisotropic behavior, 
regularity as well as stability of these models. Finally, the surface redshift 
has been calculated. All these properties of the models have been discussed 
by using the observational data of the compact stars. The plan of the paper 
is the following: In the next section, we present the review of f(T) gravity. 
Section 3 deals with anisotropic source and equation of motion equations in 
f(T) gravity with nondiagonal tetrad. Section 4 investigates with the phys¬ 
ical analysis of the proposed models. In the last section, we conclude the 
results of the paper. 

2 /(T) Gravity: Fundamentals 

In this section we briefly overview the basics of f(T) gravity (for review see ). 
The dynamical variables of TEGR are the orthonormal tetrad componenets 
&a (:r M ). The metric can be defined via these components by g, J/V = rjAB^e^, 
where t)ab — diag( 1, —1, —1, —1). Herein, A,B run over 0,1, 2, 3 for the tan¬ 
gent space of the manifold and /i, v are coordinate indices on the manifold 
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which also run over 0,1,2,3. The tetrad fields define a connection of the 
form T^ = efd^e 1 , which is named as Weitzenbock’s connection. This con¬ 
nection involves no contribution from curvature as compared to Levi-Civita 
connection for which torsion is zero. 

The torsion T* v and contorsion K^ v tensors are defined by 

T iv = e \{ d n e t (!) 

K = -\(K -T?-T?). (2) 

The torsion scalar defines the teleparallcl Lagrangian density and is given by 


X 


r-n\ qH'V 
1 flV^X ? 


(3) 


where 

ST = \{K + ( 4 ) 

One can define the modified teleparallel action by replacing X with a function 
of X, in analogy to f(R ) gravity [Tf IB], as follows 


X 


dx 4 e 


1 

2^? 


f(T) + C m , 


(5) 


where e = det(e^) = y/—g, k 2 = 87 tG = 1. The variation of the action with 
respect to the vierbein vector field ef presents 

e-'dTzST)f\T)-e\T^STl\T)+SJdTT)l'XT)+-/ A f{T) = ye'Tj"!', 

(6) 

where rj™'* is the matter energy momentum tensor and prime denotes the 
differentiation with respect to X. 


2.1 Anisotropic Matter Configuration in f(T) Gravity 

We assume the Krori and Barua spacetime to describes the strange star 
stellar configuration, which is given by [30] 

ds 2 = -e^ r) dt 2 + e u{r) dr 2 + r W, (7) 
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with v = Ar 2 , /i = Hr 2 + C where dfi 2 = dd 2 + sm 2 6d(j) 2 , A, B and C are 
arbitrary constant to be evaluated by using some physical constraints. In 
this setup we consider the off-diagonal tetrad matrix of the form j3T] 


t 

e M/ 2 0 

0 

0 

\ 


0 e v l 2 sinQcos(\) 

rcosBcoscj) 

—rsinOsincj) 



0 e v ^ 2 sin9sin(f) 

rcosdsiruj) 

rsinOcoscj) 


V 

0 e u t 2 cos6 

—rsinO 

0 

/ 


The determinant of e* can be obtained as e = e^ +l/ r 2 sinQ. 

h L 

We assume the anisotropic fluid as interior of compact object defined by 
the energy momentum tensor T^ as 

T af 3 = (p + Pt)u a up - Ptg a /3 + (Pr - Pt)v a V/ 3 , (9) 

where u a = e^<5°, v a = e%d°, p, p r and p t correspond to energy density, 
radial and transverse pressures, respectively. The dynamical equations for 
spacetime ((7j) lead to 


P 

Pr 


Pt 


T 


' "' 2 (e-' 2 - 1 )F' - ( \ - 3 - ^(1 - rv>) ) X + J ~, ( 10 ) 


F L 

2 ' 4 ' 


T 


77 - 3 + “W( 1 + r ^ / ) 77 _ T> 


2 

— 1/ 


//' 1 

77 + - 
2 r 


F f 


/ 


( 11 ) 


W2 


F' + 


T 


+ e 




( 12 ) 


where F = fx, prime denotes the derivative with respect to radial coordinate 
and T is the torsion scalar given by 


T(r) 


2e-‘'(e"'' 2 - lje"'' 2 - 1 - r/i' 
^2 


(13) 


Eqs. (ITTTjl - (TT2li represents the energy density and anisotropic pressures for the 
off-diagonal tetrad. In case of diagonal tetrad there exist an off-diagonal 
equation which results in linear form of algebraic function f(T). However, 
the choice of off-diagonal tetrad does not provide any particular condition 
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and one can set a consistent f(T) model. In this setup we consider a obser- 
vationally viable f(T) model given by [32, 33] 

f(R) = PT n , (14) 


where (5 is an arbitrary constant. In [33] Banrba et al. discussed the existence 
of finite-time future singularities for power law type model of the form T 13 , 
it is found that the terms involving powers like (3 > 1 can remove such 
singularities. Setare and Darabi [34] developed the f(T) model for power law 
solutions and discussed the case when universe enters the phantom phase. 

Using the f(T) model (1I2|) together with relations of p and v in metric 
(J4J) , we find the expressions for p, P r and Pt in the following form 


P = 


-)Tl —2 


P 


,Ar 2 


2\2 


(1 + 2 Br 2 ) 


( e Ar*/2 _ !)( e ArV2 _ 2 Br 2 - l) 2 

4n 2 (l + Ar 2 + 2 ABr 4 ) + 2n{2 + 3 Ar 2 + Br 2 + 6 ABr 4 + 2 B 2 r A ) 


+ e Ar (-3 - An 2 - ABr 2 + 2n(3 + Br 2 )) + e Ar /2 (3 + 8 Br 2 + 25 2 r 4 
+ 4n 2 (2 + Ar 2 + ABr 4 ) — 2n(5 + 3 Ar 2 + 4i?r 2 + 2ABr 4 + 2 B 2 r 4 ))) 


x 


; _ A , 2 ( e Ar 2 /2 _ !)( e Ar 2 /2 _ 2Br 2 - 1)' 


(15) 


'in—2 o„—Ar 2 


Pr = 


2 pe~ 


; -Ar 2 ( e Ar 2 /2 _ l^Ar 2 /2 _ 2 Br 2 - 1) ' 


71—1 


^Ar 2 


+ 2e Ar2/2 (n - 1)(1 + Br 2 - (2 n - 1)(1 + 2 Br 2 ))} , 


(16) 


Pt = r~ 4 /3 


-A + 2n 2 


(n — 1 )Ar 2 (n — 1 )r 2 (A — AB + 2ABr 2 ) 


(e Ar2 /2 - 1)2 ( e Ar2/2 _ 2Br 2 - l) 2 

An — 2 B + B(An — B)r 2 2£> + (A — 4£>)n + B(An — B)r 2 


x 


B(e Ar2 / 2 - 1) 

l)( e Ar 2 /2 


3 —Ar 2 ( p Ar 2 /2 


B(e Ar2 / 2 
2Br 2 - 1)' 


2 Br 2 - l) 2 


(17) 


Using the expressions of energy density, radial and transverse pressures (1T5]) - 
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m, we find the corresponding equation of state (EoS) as follows 


(1 - e Ar2/2 + 2 Br 2 ){e Ar2 + 2 e Ar2/2 {n - 1)(1 + Br 2 ) - (2 n - 1)(1 


+ 

+ 

+ 

+ 


UJ t = 


X 

+ 

+ 

+ 

+ 


3r 2 ))] / e 3Ar2/2 - (1 + 2 Br 2 ) 2 - 4n 2 (l + Ar 2 + 2 ABr 4 ) + 2n{2 

1 r 2 + Br 2 + 6 ABr 4 + 2B 2 r 4 ) + e Ar2 (-3 - 4n 2 - 4 Br 2 + 2n(3 
r 2 )) + e Ar2/2 (3 + 8 Br 2 + 45 2 r 4 + 4n 2 (2 + Ar 2 + A5r 4 ) - 2n(5 
lr 2 + ABr 2 + 2ABr 4 + 25 2 r 4 ))] , (18) 

2Ar2 (n - 1) + 2 e 3Ar2/2 (n - l)(n - 2 - 25r 2 ) - (1 + 25r 2 ) 2 - 2n 2 

+ 5r 2 )(l + A(r 2 + 2 -Br 4 )) + n(3 + 85r 2 + 9 5 2 r 4 + 25 3 r 6 + A(r 2 
3r 4 + 25 2 r 6 )) + e^ r2 (-2(3 + 65r 2 + 25 2 r 4 ) - 2n 2 (3 + Br 2 + A(r 2 
r 4 ))+n(12 + 165r 2 + 55 2 - 4 w ° ~ 

3 2 r 4 — n( 5 + 105r 2 + 75' 

- 2 f9 4- AUr 2 4- R 2 'T 4 illl / 


3 + 65r + 25 r ) - 2rA(3 + 55 + A(5 
,2 r 4 + A(r 2 + Br 4 ))) + 2e Ar2/2 {2 + 65r 2 
! r 4 + 5 3 r 6 + A(r 2 + 2 5r 4 )) + n 2 (3 + 25r 2 


+ Ar 2 + 2A5r 4 ) + 2n(2 + 3 Ar 2 + 5r 2 + 6A5r 4 + 25 2 r 4 ) + e Ar2 (-3 
— 4n 2 — 45r 2 + 2n(3 + Br 2 )) + e Ar2 ^ 2 (3 + 8 Br 2 + 45 2 r 4 + 4n 2 (2 + Ar 2 
+ ABr 4 ) — 2n(5 + 3 Ar 2 + 45r 2 + 2 ABr 4 + 25 2 r 4 ))] . (19) 


3 Physical Analysis 

In this section we explore some features of the anisotropic compact star which 
includes anisotropic behavior, matching conditions and stability constraints. 

3.1 Anisotropic Behavior 

Here, we present the evolution of p, p r and pt as shown in Figs. 1-3 for the 
strange star candidates Her X-l, SAX J 1808.4-3658 and 4U 1820-30 (see 
table 1). Herein, we set n = 2 and /3 = —2. From Eqs. (1T5T) and (115]) . 
one can End the expressions for dp/dr, d?p/dr 2 , dp r /dr and d 2 p/dr 2 . We 
show the results for dp/dr and dp r /dr in Appendix as Eqs. flA. II) and (]A.2|) 
respectively. Figures 4 and 5 show the variation of derivatives of p and p r 
with respect to radial coordinate for strange star Her XI. We hnd that 
dp/dr < 0 and dp r /dr shows transition from positive to negative values. 
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Figure 1: Evolution of energy density p versus radial coordinate r(km ) at 
the stellar interior of strange star candidates. 



Figure 2: Evolution of radial pressure p r versus radial coordinate r(km) at 
the stellar interior of strange star candidates. 
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Figure 3: Evolution of tangential pressure p t versus radial coordinate r(km ) 
at the stellar interior of strange star candidates. 
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Figure 4: Evolution of dp/dr and d 2 p/dr 2 versus r(km ) at the stellar interior 
of Her XI. 
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Figure 5: Evolution of dp r /dr and d 2 p r /dr 2 versus r(km ) at the stellar inte¬ 
rior of Her XI. 
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Figure 6: Evolution of ay versus radial coordinate r(km ) at the stellar interior 
of strange star candidates. 



Figure 7: Evolution of radial pressure c o t versus radial coordinate r{km ) at 
the stellar interior of strange star candidates. 


We also examine behavior of derivatives at r = 0, centre of strange star 
ffer XI. It can be seen that 


dp dp r 

dr ’ dr 


d 2 p 

dr 2 


0, 


d 2 p r 

dr 2 


> 0 . 


( 20 ) 


Clearly, at r = 0, density is maximum and it decreases outward with the 
increase in r as shown in Figure 1. In case of p r , we have increasing func¬ 
tion unlike the evolution p, which attains maximum value in midway to the 
boundary of star and then decreases. One can test the evolution of p and p r 
for the strange stars SAX J 1808.4-3658 and 4U 1820-30, which would show 
similar behavior. The evolution of EoS parameters ay and ay is also presented 
for the strange stars in f(T) gravity. It can be seen that ay > 0, which fa¬ 
vors the quintessence regime near the boundary of star whereas 0 < ay < 1 
despite of the role of modified gravity. 
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Figure 8: Evolution of anisotropy measurement A versus radial coordinate 
r(km ) at the stellar interior of strange star candidates. 


The anisotropy measurement in f(T) gravity, for spherically symmetric 
metric is given by 


A = 


x 


2 (Pt-Pr) ri2 n f _ Ar 2 // Arl 


r 

2 Ar 2 


2tx—|— 1 


-1 


- 1 - 2 Br‘ 


n —2 


3 A- 


+ 2e 2 (—2 + n — Br 2 ) + e Ar ' 2 (6 + Ar 2 + 4£>r 2 + ABr 


+B 2 r 4 — 2n(3 + Ar 2 + Br 2 + ABr 4 )) + (1 + Br 2 )( 1 + Ar 2 — Br 2 

+2AE>r 4 + 2 B 2 r 4 — 2n(l + Ar 2 + 2 ABr 4 )) + 2e~ (—2 — Ar 2 — Br 2 
-2ABr 4 - B 2 r 4 - B 3 r 6 + n(3 + 2 Br 2 + Ar 2 (2 + ABr 2 + B 2 r 4 )))} (21) 


The measure of anisotropy depends on radial and tangential pressures, it is 
directed outward when p t > p r which implies A > 0 whereas p t < p r results 
in A < 0 so that anisotropy is directed inward. We present the evolution of 
A for strange star candidates Her X-l, SAX J 1808.4-3658 and 4U 1820-30 in 
Fig. 8. It can be seen that A > 0 for our model which implies the existence 
of repulsive force, allowing more massive distribution. One can see that f(T) 
gravity does not affect the anisotropic force. Here, A vanishes at the center 
r = 0. 


3.2 Matching Conditions 

Here, we match the interior metric ([7j) to the vacuum exterior spherically 
symmetric metric given by 

ds 2 = — ^1 — —— ^ dt 2 + ^1 — —— ^ dr 2 + r 2 d6 2 + r 2 Sin 2 6dip 2 , (22) 
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At the boundary surface r = R continuity of the metric functions g u , g r r and 
yield, 


9tt 9tt ? 


9rr 



dg tt _ dgit 

dr dr ’ 


(23) 


where — and +, correspond to interior and exterior solutions. From the 
interior and exterior metrics, we get 


A 

B 

C 




1 - 

2 M 
~R 
2 AT 



2M\"' 

~r) 


(24) 

(25) 

(26) 


For the given values of M and R for given star, the constants A and B takes 
the values as given in the table 1. 


Table 1: Values of constants for given Masses and Radii of Stars [351136| 


Strange Quark Star 

M 

R(km ) 

M 

R 

A(km 2 ) 

B(km 2 ) 

Her X-l 

0.88 M e 

7.7 

0.168 

0.006906276428 

0.004267364618 

SAX J 1808.4-3658 

1.435Af 0 

7.07 

0.299 

0.01823156974 

0.01488011569 

4U 1820-30 

2.25 M q 

10.0 

0.332 

0.01090644119 

0.009880952381 


3.3 Energy Conditions 

Energy conditions have been a handy tool to limit the arbitrariness in the 
energy-momentum tensor, based on Raychaudhuri equation with attractive¬ 
ness property of gravity. These conditions include weak energy condition 
(WEC), null energy condition (NEC), strong energy condition (SEC) and 
dominant energy condition (DEC). In case of anisotropic fluid, we find the 
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Figure 9: Evolution of energy constraints at the stellar interior of strange 
star Her X-l. 


following inequalities 


NEC : 

P + Pr 

WEC : 

P > o, 

SEC : 

P + Pr 

DEC : 

P > \P 


> 0, p + p t > 0, 

P + Pr > 0, p + 

> 0, p + p t > 0, 

■I, P>\Pt\- 


Pt > 0, 

P + p r + 2p t > 0, 


We examine these constraints for the strange star Her X-l and show the 
evolution in Figure 9. ft can be seen that these conditions are satisfied for 
strange star Her X-l and one can also the constraints on other strange star 
models. 


3.4 TOV Equation 

For an anisotropic fluid the generalized Tolman-Oppenhcimer-Volkoff (TOV) 
equation has the following form 

dPr + v'(p + Pr) + 2(p r ~ Pt) = Q 
dr 2 r 

Following [37j, we can express the above equation in terms of gravitational 
mass and henceforth it results in the equilibrium condition for the strange 
star, involving the gravitational, hydrostatic and anisotropic forces of the 
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Figure 10: Evolution of gravitating, hydrostatic and pressure anisotropic 
forces at the stellar interior of strange star Her X-l. 

stellar object as 


F g + Fh + F a — 0, 

F g = -Br{p + p r ), F h = -^, F a 


2(Pt ~ Pr) 
r 


(28) 


Using the effective p, p r and p t (TTOh - (TT2T) . we can find the numerically the 
equilibrium condition in f(T) theory. In Figure 10, we show the evolution 
of above forces at the interior of strange star Her X-l. This figure indicates 
that the static equilibrium can be attained due to pressure anisotropy, grav¬ 
itational and hydrostatic forces for strange star candidates in f(T) gravity. 
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3.5 Stability Analysis 

In this section, we discuss the stability of strange star models in f(T) theory. 
The radial and transverse speeds can be calculated as 


v 


2 

sr 


(-1 + e 6 ^ - 2 Br 2 )(e 2Ar2 + 2e^(ra - 2)(1 + Br 2 ) - (2 n - 1)(1 


+2Br 2 )(l + A(r 2 + 2 Br 4 )) + eUr 2 )( 6 + 10Br 2 + A(r 2 + 2 Br 4 + 2 B 2 r 6 ) 
—2n(3 + 4 Br 2 + A(r + Br 3 ) 2 )) + 2e^(-(l + 2Br 2 )(2 + A(r 2 + .Br 4 )) 


+n(3 + 5 Br 2 + 2A(3Br 4 + r 2 


+ 2i? 2 r 6 ))))]/ 


5Ar 2 

1 — e 2 


+ 7Ar 2 — 10 Br 2 


+18 ABr 4 + 8 B 2 r 4 + 12AB 2 t 6 + 8 AB 3 r 8 + e 2Ar \h + 4n 2 - 2Ar 2 + 6 Br 2 
+2n (—3 + Ar 2 — Br 2 )) + 4n 2 (l + A(r 2 + 2 Br 4 )) 2 — 2n(2 — 7 Br 2 + 2B 2 r 4 

+3 r 4 (A + 2ABr 2 ) 2 + Ar 2 (7 + 15.Br 2 + 45 2 r 4 + 45 3 r 6 )) - e^(10 + 4 Br 2 
+8 B 2 r 4 - 2A 2 r 4 (3 + 2 Br 2 ) + Ar 2 (16Br 2 + 2B 2 r 4 + 3) + 8n 2 (2 + A(r 2 
+Br 4 )) + 2n (—11 + 2Br 2 — 2B 2 r 4 — Ar 2 ( 5 + 12 fir 2 + £> 2 r 4 ) + 2A 2 (r 4 
+5r 6 ))) + e Ar2 (-4^V(3 + 55-r 2 ) + 2(5 - 10Br 2 + 8 B 2 r 4 ) + Ar 2 (19 
+56B>r 2 + 16B 2 r 4 + 4£> 3 r 6 ) + 4n 2 (6 + r 4 (A + ABr 2 ) 2 + A(6r 2 + 8 Br 4 )) 

— 2n(15 — 14.Br 2 + QB 2 r 4 + hl 2 r 4 (—1 — 2Br 2 + 2 B 2 r 4 ) + 2Ar 2 (21Br 2 

+4 B 2 r 4 + B 3 r 6 + 10))) + e^(-5 + 28 Br 2 - 165 2 r 4 + 27L 2 r 4 (3 + 8 Br 2 
+4 B 2 r 4 ) — Ar 2 (21 + 58 Br 2 + 30 B 2 r 4 + 12B : V 6 ) — 8n 2 (2 + Ar 2 ( 3 + 5 Br 2 ) 
+A 2 (r 4 + 3 Br 6 + 2 B 2 r 8 )) + 2n(9 — 18.Br 2 + 6 B 2 r 4 + Ar 2 (21 + 45£>r 2 
+135 2 r 4 + 6B 3 r 6 ) + 4A 2 (r 4 + 35r 6 + 25 2 r 8 )))] . (29) 
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v 2 t — — (—1 + e 3 2 ’ (n — 1) — 2 Ar 2 — 10 ABr 4 — 6 B 2 r 4 — 18A5V — 65 3 r 6 

—12 A5 3 r 8 — 4 5 4 r 8 + e^~ (n — 1)(—6 + 2 n — 5r 2 ( 4 + Ar 2 )) — 2n 2 (l 
+5r 2 )(l + A(r 2 + 2 5r 4 )) 2 + n(3 + 25r 2 + 5 5 V + 2B 3 r 6 + r 4 (l + 5r 2 ) 
x(A + 2ABr 2 ) 2 + Ar 2 (6 + 25 Br 2 + 395V + 28 B 3 r 6 + 45 4 r 8 )) - e 2Ar2 (15 
+165r 2 + 6 5 2 r 4 + Ar 2 (2 + 105r 2 + 3 B 2 r 4 ) — A 2 (r 4 + Br 6 ) + 2n 2 (5 + Br 2 
+2A(r 2 + Br 4 )) + n (—25 — 185r 2 — 5 B 2 r 4 — Ar 2 (6 + 17 Br 2 + 45 2 r 4 ) + 2A 2 (r 4 

+Br 6 ))) + e^(20 + 245r 2 + 24B 2 r 4 + 6B 3 r 6 - 3A 2 (r 4 + 2 J Br 6 ) + Ar 2 (UBr 2 
+8 + 15i? 2 r 4 + B 3 r 6 ) + 2n 2 (10 + 45r 2 + r 4 (A + ABr 2 ) 2 + 2Ar 2 (6Br 2 + B 2 r 4 
+4)) + n(A 2 r 4 ( 5 + 10i?r 2 + 2B 2 r 4 ) — 2(20 + 16£>r 2 + 10i? 2 r 4 + B 3 r 6 ) — Ar 2 (24 
+70fir 2 + 355 2 r 4 + 5 B 3 r 6 ))) - e Ar2 (15 + 165r 2 + 36 B 2 r 4 + 185 3 r 6 + 45 4 r 8 
-3A 2 (r 4 + 3 Br 6 + 25V) + Ar 2 (12 + 52 Br 2 + 395 2 r 4 + 65 3 r 6 - 25 4 r 8 ) 

+2?r 2 (10 + 6 5r 2 + 4Ar 2 (3 + 65r 2 + 25 2 r 4 ) + A 2 r 4 (3 + 9 5r 2 + 75 2 r 4 + 5 3 r 6 )) + 
n(—35 — 28 Br 2 — 30 5 2 r 4 — 65 3 r 6 + A 2 r 4 (3 + 9 Br 2 + 85 2 r 4 + 25 3 r 6 ) — Ar 2 (36 

+1185r 2 + 975 2 r 4 + 305 3 r 6 + 25 4 r 8 ))) + e^(6 + 45r 2 + 245V + 185 3 r 6 
+85 4 r 8 + Ar 2 (8 + 375r 2 + 455 2 r 4 + 175 3 r 6 - 25 4 r 8 ) - A 2 (r 4 + 45r 6 + 65 2 r 8 
+45 3 r 4 0) + 2n 2 (5 + 45r 2 + 2Ar 2 (4 + 105r 2 + 55 2 r 4 ) + AV(125r 2 + 145 2 r 4 
+3 + 45 3 r 6 )) + n(—2(8 + 65r 2 + 105 2 r 4 + 35 3 r 6 ) + A 2 r 4 (—1 — 45r 2 — 25 V 

+45V) - Ar 2 (24 + 895r 2 + 1055V + 535 3 r 6 + 65 4 r 8 ))))/((-l + e^(l 

+ 7Ar 2 - 105r 2 + 18A5r 4 + 85 2 r 4 + 12A5V + 8A5 3 r 8 + e 2Ar2 (5 + 4n 2 
—2Ar 2 + 65r 2 + 2?r(—3 + Ar 2 — Br 2 )) + 4n 2 (l + A(r 2 + 25r 4 )) 2 — 2n(2 — 7 Br 2 

+25 2 r 4 + 3r 4 (A + 2ABr 2 ) 2 + Ar 2 (7 + 15 Br 2 + 45 2 r 4 + 45 3 r 6 )) - e^(10 
+45r 2 + 85 2 r 4 — 2A 2 r 4 (3 + 2 5r 2 ) + Ar 2 (3 + 165r 2 + 25 2 r 4 ) + 8n 2 (2 + A(r 2 
+5r 4 )) + 2n(—11 + 2 Br 2 - 25V - Ar 2 (5 + 125r 2 + 5 2 r 4 ) + 2A 2 (r 4 + Br 6 ))) 
+e Ar2 (—4A 2 r 4 (3 + 5 Br 2 ) + 2(5 — 105r 2 + 85 2 r 4 ) + Ar 2 (19 + 56 Br 2 + 165 V 
+45 3 r 6 ) + 4n 2 (6 + r 4 (A + ABr 2 ) 2 + A(6r 2 + 8 Br 4 )) — 2n(15 — 145r 2 + 65 2 r 4 

+AV(-1 - 25r 2 + 25V) + 2Ar 2 (10 + 215r 2 + 45V + 5 3 r 6 ))) + e^(-5 
+285r 2 — 165 2 r 4 + 2AV(3 + 8 Br 2 + 45V) — Ar 2 (21 + 58 Br 2 + 305 2 r 4 
+125 3 r 6 ) — 8n 2 (2 + Ar 2 (3 + 5 Br 2 ) + A 2 (r 4 + 3 Br 6 + 2 B 2 r 8 )) + 2n(9 — 185r 2 
+65 2 r 4 + Ar 2 (21 + 45 Br 2 + 135V + 65 3 r 6 ) + 4A 2 (r 4 + 3 Br 6 + 25 2 r 8 ))))). (30) 
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Figure 11: Evolution of v 2 r versus radial coordinate r(km ) at the stellar 
interior of strange star candidates. 



Figure 12: Evolution of v 2 t versus radial coordinate r{km ) at the stellar 
interior of strange star candidates. 


In [38J, Herrera developed a new technique to explore the potentially unsta¬ 
ble matter configuration and introduced the concept of cracking. One can 
analyze the potentially stable and unstable regions regions depending on the 
difference of sound speeds, the region for which radial sound speed is greater 
than the transverse sound speed is said to be potentially stable. In Figure 
11 and 12, we plot the radial and transverse speeds for different strange star 
candidates. Here, the values of v 2 sr and v 2 t lies in the range 0 < \v 2 \ < 1 
within the anisotropic matter configuration. We also present the evolution 
of v 2 t — v 2 r in Figure 13, which shows that difference of two sound speeds 
satisfies the inequality \v 2 t — v 2 r \ < 1. Moreover, v 2 t — v 2 r retain the similar 
sign for Her X-l whereas the sign changes in case of SAX J 1808.4-3658 and 
4U 1820-30. Hence, our proposed strange star model is stable for Her X-l 
and unstable for both SAX J 1808.4-3658 and 4U 1820-30. 
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Figure 13: Evolution of v 2 t — v 2 r versus radial coordinate r(km ) at the stellar 
interior of strange star candidates. 

3.6 Surface Redshift 

The compactness of star in f(T) gravity can be found as 
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The surface redshift corresponding to compactness is given by 
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Figure 14: Evolution reshift Z s versus radial coordinate r(km ) at the stellar 
interior of strange star Her X-l. 


In Figure 13, we show the redshift plot for the star Her X-l and maximum 
redshift turns out to be Z q = 0.025. 


4 Conclusion 

In this paper, we have addressed the issue of modeling the spherically anisotropic 
compact stars in f(T) gravity using the off diagonal tetrad field. For the an¬ 
alytic solution of the system of differential equations, we have assumed the 
KB form of the metyric components and the power law form of f(T) model, 
as f(T) = (3T n , with n ^ 1 and /3 is an arbitrary small constant, which plays 
the role of departure of these theories from GR. This approach leads to the 
set field equations, which second order nonlinear differential equations, for 
the analytic form of solutions we assume that the metric coefficients for the 
interior and exterior regions of a star continuous over the boundary. This 
is the smooth matching of the general spherical interior geometry with the 
Schwarzschild exterior geometry, because in both regions metric function are 
continuous upto fist order of derivatives (GR and f(T) are second order 
derivative theories). This matching helps to use the observational data of 
stars for the calculation of arbitrary constants appearing in the system due 
to parameterization of metric functions. 

We have determined the explicit form of the matter density and pressure 
components, which are regular inside the stars as shown in figures (1-3). The 
behavior of EoS parameter implies that near the boundary the matter inside 
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the stars behaves as quintessence in the presence of f(T) terms (see figure 6). 
The anisotropic parameter A = 0 at center vanished and it becomes A > 0, 
everywhere inside the star except at center, which predicts that there may 
exist the repulsive force which leads to the formation of more massive star. 
The physical viability of the solutions lies in the fact that solutions satisfies 
the energy conditions (Figure 9). The analysis of TOV equation presented 
in this paper shows that the static equilibrium can be achieved for strange 
stars in f(T) as shown in figure 10. The stability of the proposed model is 
given in figures 11-13, on the physical investigations of these figures predict 
that our proposed strange star model is stable for Her X-l, and unstable for 
both SAX J 1808.4-3658 and 4U 1820-30. The maximum redshift of Her X-l 
turn out be 0.025 which is shown in Figure 14. 
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x(10 — 20 Br 2 + 16B 2 r 4 — 4A 2 r 4 (3 + 5 Br 2 ) + Ar 2 (56Br 2 + 16B 2 r 4 + 4B 3 r 6 
+19) + 4n 2 (6 + r 4 (A + ABr 2 ) 2 + A(6r 2 + 8 Br 4 )) — 2n(15 — 14Br 2 + 6 B 2 r 4 

+A 2 r 4 (—1 — 2 Br 2 + 2 B 2 r 4 ) + 2Ar 2 (10 + 21 Br 2 + 4B 2 r 4 + B 3 r 6 ))) + e~ (—5 
+28Br 2 — 16B 2 r 4 + 2A 2 r 4 (3 + 8 Br 2 + 4B 2 r 4 ) — Ar 2 (58Br 2 + 30B 2 r 4 + 12B 3 r 6 
+21) - 8n 2 (2 + Ar 2 (3 + 5Br 2 ) + A 2 (r 4 + 3Br 6 + 2B 2 r 8 )) + 2n(-18Br 2 + 9+ 
6B 2 r 4 + Ar 2 (21 + 45Br 2 + 13B 2 r 4 + 6B 3 r 6 ) + 4A 2 (r 4 + 3Br 6 + 2B 2 r 8 )))] ,(A.l) 
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+1) + 2n 2 (5 + 4Sr 2 + 2Ar 2 (4 + lOSr 2 + 5S 2 r 4 ) + A 2 r 4 (12Sr 2 + 14S 2 r 4 + 4S 3 r 6 
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